Quantum phase transitions in the Fermi-Bose Hubbard model 
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We propose a multi-band Fermi-Bose Hubbard model with on-site fermion-boson conversion and 
general filling factor in three dimensions. Such a Hamiltonian models an atomic Fermi gas trapped 
in a lattice potential and subject to a Feshbach resonance. We solve this model in the two state 
approximation for paired fermions at zero temperature. The problem then maps onto a coupled 
Heisenberg spin model. In the limit of large positive and negative detuning, the quantum phase 
transitions in the Bose Hubbard and Paired-Fermi Hubbard models are correctly reproduced. Near 
resonance, the Mott states are given by a superposition of the paired-fermion and boson fields and 
the Mott-superfluid borders go through an avoided crossing in the phase diagram. 

PACS numbers: 



The experimental investigation of cold atomic gases 
is proceeding rapidly. Both bosons and fermions have 
been brought to quantum degeneracy P, . They have 
been trapped in the sinusoidal lattice potential created 
by an optical standing wave of two counter-propagating 
lasers 0, 0]. In the tightly bound regime of this po- 
tential, the Bose-Hubbard Hamiltonian has proven to 
be a useful model to describe the transition from a su- 
perfluid, in which the atoms are delocalized, to a Mott 
insulator, which has an integer number of atoms at 
each lattice site |M S 0- Recently, the success- 
ful implementation of Feshbach resonances in degenerate 
Fermi gases has enabled the experimental study of the 
Bardecn-Cooper-Schrieffer (BCS) to Bose-Einstein con- 
densate (BEC) crossover in the continuum, a long stand- 
ing theoretical problem 0, 0| . Initial evidence of a new 
superfluid state has been found in the strongly interact- 
ing regime 0. It is a logical next step to study such 
a crossover in an atomic Fermi gas trapped in a lattice 
potential 0. 

In this Letter, we investigate the BCS-BEC crossover 
in the context of the Fermi-Bose Hubbard Hamiltonian 
(FBHH), motivated by this vigorous experimental activ- 
ity. Hubbard models have proven useful in experiments 
on BEC's [^11] and are expected to be equally relevant 
for fermions 0, ll^ . A phenomenological fermion-boson 
conversion term in a simplified FBHH was first suggested 
in the context of high temperature superconductors , 
while FBHH's without conversion have been treated in 
the context of cold quantum gases 0|. Unlike in this 
earlier work, the model we shall study includes the possi- 
bility of both classical and quantum phase transitions Q 
in the Fermi and Bose Hubbard limits, as well as a con- 
version term. In addition, we allow the fermions to oc- 
cupy multiple bands, so that the filling factor is not con- 
strained. In contrast to high- Tc ^^li fermion-boson con- 
version is a real physical process in cold quantum gases, 
where a Feshbach resonance is used to coherently transfer 
fermionic atoms into a bound two-atom bosonic state, as 
illustrated in Fig. ^ 



The effect of the conversion term is to lock the order 
parameter of the fermions and bosons together. It thus 
leads to a reduction in the number of quantum phases 
from four to two. The main reason we introduce a bosonic 
field is to describe the the BCS-BEC crossover: the at- 
tractive Fermi-Hubbard Hamiltonian, even in the paired 
fermion limit, does not map simply onto the repulsive 
Bose Hubbard Hamiltonian After proposing this 

new FBHH, we solve it in detail in the limit of on-site 
paired fermions in the two-state approximation at 
zero temperature for a filling of from zero to two fermions 
per site. The on-site paired-fermion limit corresponds to 
the experimentally realizable case of a strongly confining 
potential and/or strong interactions. In this limit the 
FBHH maps isomorphically onto a coupled Heisenberg 
spin model, or coupled magnets. 

Consider the FBHH in the grand canonical ensemble. 
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Equations (0)-® are the usual repulsive Bose-Hubbard 
and multi-band attractive Fermi-Hubbard Hamiltonians 
for a uniform lattice and Eq. is the fermion-boson cou- 
pling. The symbol (i, j) denotes nearest neighbors, while 
the indices s G {Ti i} ^^nd m denote the spin state and 
band number. The hopping or tunnelling strengths J/^f, 
and the on-site interaction strengths Vf.b are taken as real 
and positive definite. The band gap energy of the to'^ 
band is Em ■ The strength g of the interconversion term 
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FIG. 1: Outer figure: pairs of fermionic atoms in the 
open channel are coherently transferred into a closed chan- 
nel, bosonic state via a Feshbach resonance. Inset; second 
order degenerate perturbation theory in the limit Jf <^ Vj 
leads to two hopping events on the lattice, (a) pair hopping, 
and (b) a single fermion hopping to an adjacent site and back. 



and Vfb of the density coupling may have either sign. 
The creation and annihilation operators /^,/ and b\b 
satisfy the usual commutation relations for fermions and 
bosons, respectively. The number operators are defined 
as n'^i = bjbi, n{^^ = f-^^fism- In order to match the 
physical context of quantum degenerate gases in chemi- 
cal equilibrium, we require 

fib = 2^/ + fiv , (5) 

where u is the detuning associated with a Feshbach res- 
onance and we set h — 1. The conserved quantity 

^^2E,".? + E.,.,™"L„ (6) 

is the total number of fermions. Eliminating fib by sub- 
stituting Eq. 10 into Eqs. Q-Q), one finds that fif mul- 
tiplies n. One can thus take fif as the chemical potential 
of the coupled system, while v determines the relative 
number of bosons and fermions. 

The FBHH of Eqs. models a pscudo-spin-1/2 

system of fermions with s-wave interactions, as in exper- 
iments IMSOl- In practice, the index s G {T7i} rep- 
resents two hypcrfine states in the level structure of an 
effectively fermionic alkali atom, such as ''"K or ^Li, scat- 
tering near threshold in an open channel. The bosonic 
field represents a bound closed-channel molecular state, 
^Li2 or ^°K2, which is coupled to the fermionic field via 
a resonance with an unbound open-channel atomic state, 
called a Feshbach resonance. A schematic is shown in 
Fig. n^a). Note that Vf^b and g are not functions of 
1^. Methods for calculating the parameters Vf,Vb, etc. 
in Eqs. from few-body atomic physics have been 

described in detail elsewhere |l7l| . Another important 



assumption is that the pairing of fermions into bosons 
occurs on-site. This is physically reasonable for present 
experiments (lif . 

We consider the limit in which Jf <^ Vr, which cor- 
responds to a strongly confining lattice Since the 
lattice height is proportional to the intensity of the lasers 
creating the standing wave, this is straightforward to 
obtain. Because the on-site interactions are attractive 
and s-wavc, and the hopping is taken perturbatively, 
the fermions form spin- up/spin-down pairs. We also re- 
strict them to be in the lowest band. This is the typical 
experiment case in three dimensions, where 10^ to 10^ 
fermions are distributed among 100^ sites. Thus m = 1 
and n G [0, 2], i.e., there are from zero to two fermions, or 
zero to one fermi pair, per site. Then second order degen- 
erate perturbation theory maps Hf onto a new spin-1/2 
system (a quantum XXZ model) [l5| : 

where Jj- = 8Jj/Vf and fi'j = 2fif + Vf/2. The operator 
= (T~)t = //|/,/| is a pair creation/annihilation oper- 
ator and Ui = + — 1). The perturbative treat- 
ment results in two hopping-type events, as is sketched 
in the inset of Fig. the T^^'r," term corresponds to 
pair hopping, while the niHj term corresponds to a single 
fermion hopping to an adjacent site and hopping back. 

These operators obey the commutation relation 
[rf,r/] = 2iT^e"l'^6,j, where a,/?, 7 € {x,y,z}, r± = 
rf ± irf , and = n^. Thus, despite the fact that is 
a creation/annihilation operator for fermion pairs, the t 
operators obey the Pauli spin commutation relations, not 
the bosonic commutation relations. This is one reason 
why the attractive Fermi-Hubbard model does not map 
simply onto the repulsive Bose-Hubbard model, even in 
the limit of strong interactions. A second reason is that in 
order to achieve such a mapping, a sum over many bands 
is required, since the internal energy of bosons composed 
of two fermions is much greater than the band spacing. 
In contrast, the FBHH is asymptotically able to represent 
both the attractive Fermi-Hubbard and repulsive Bose- 
Hubbard models in a simple way. It is therefore a good 
candidate for the study of the BCS-BEC crossover. 

In general, a paired Fermi Hubbard Hamiltonian can 
act on all number states of the fermions. However, as 
in Eq. we consider only n S [0,2], the Hilbert space 
on which it operates is restricted to two paired-number 
states. Thus H'j: is equivalent to the Hcisenbcrg spin 
Hamiltonian, or a magnet, 

where fi'f plays the role of the magnetic field /i^. One 
therefore expects paramagnetic and either fcrro- or anti- 
ferromagnetic phases. The former correspond to the su- 
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perfluid phase, while the latter are Mott and charge- 
density wave (checkerboard) phases. Similarly, the re- 
striction of the Hilbert space on which Hh operates to 
two number states leads to an isotropic Heisenberg spin 
Hamiltonian (the quantum XX model We formu- 

late the two-state approximation for the coupled model 
as superposition states of the form It/j) = Y[j H')j, where 

|V^)j = |0)5 ® |0)^^ cos% + sin(0,) e"^^ 
X (|1)5 ® |0)J cosxj + |0)5 ® sinx, e'^"^ ) . (9) 

The superscripts b and / refer to Fock states of bosons 
and fermi pairs on the j*^ site. 

The two state approximation is useful in determining 
the Mott-superfluid borders in the phase diagram. The 
Mott state is a single number state, while the lowest order 
approximation of a superfluid is a superposition of two 
number states. Therefore, Mott states occur in Eq. Q 
for 9 £ {0,7r/2,7r}. The mixing angle x is determined 
by the detuning ly in Eq. (jSJ. To determine which phase 
is energetically favorable one evaluates £^gs = {ip\H\ip). 
We make the uniform approximation 9j = 9 = (f), 
Xj = X I ckj ~ Q:- Then (/) does not appear in the ground 
state energy, while a can only change the sign of g. Set- 
ting g' = min[(7 exp(ia)], neither nor a need be consid- 
ered to obtain the phase diagram. An important point 
is that the ground state is either paramagnetic (super- 
fluid) or ferromagnetic (Mott). It can be proven that it 
is not antiferromagnetic (charge-density wave) , either by 
setting the angles to differ by 7r/2 on each site, or by 
making a spin rotation in the Hamiltonian |20l| . 

The Mott-superfluid borders are obtained as follows. 
The ground state energy is expanded around the Mott 
angles 9 G {0, 7r/2,7r}. The zeroth order term gives the 
energy. The first order term is zero, showing that the 
Mott state is always an extremum. The sign of the sec- 
ond order term determines whether the Mott state is a 
maximum or a minimum. Setting this equal to zero, one 
obtains the Mott-superfluid borders. One must also ex- 
tremize in the mixing angle x and determine whether or 
not it is a maximum. Thus there are three conditions: 

d^Eg^/d9^ = 0, (10) 
dE^^/dx = 0, (11) 
d^E,s/dx' > 0. (12) 

Using conditions H10() - (|ll|l to eliminate x and Eq. lO to 
eliminate fib, one finds a quartic equation in /if. The 
coefficients are functions of J'j, Jb, Vj, Vt, \g'\ = \g\, 
and h'. The solution to this quartic equation, though 
lengthy, can be written in closed analytic form. It is best 
understood when evaluated in limits of the parameters 
and for particular values of them. 

First consider the case v — > ±oo. We assume a bi- 
partide lattice with Z the number of nearest neighbors. 
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FIG. 2: (color online) Shown is the phase diagram for de- 
tunings (a) u/Vf = -10, (b) u/Vf = -1, (c) u/Vf = 1/2, 
(d) ly/Vf = 10. The blue solid curves show the Mott- 
superfluid borders, while the red dashed curves show alter- 
nate extrema which are maxima (see Fig.|^. SF = superfluid, 
MI = Mott insulator, n = fermion filling factor. 



Then x G {0, 7r/2} and one obtains the Mott borders 

- -l/4+(Z/2)(l-2a/)4/I/;, (13) 
fib/Vb - -2abZJb/Vb, (14) 

where ct/ = ±1 gives the vacuum/one- fermi-pair and 
(Tf, = ±1 the vacuum/one-boson Mott states. Equa- 
tions (|13|I - H14(I correspond to the solutions one finds for 
g = in the two-state approximation. For ly — > -|-oo, con- 
dition (c) shows that Eq. H13(l is a minimum and Eq. H14|) 
is a maximum. For ly — oo, the inverse is the case. 
Thus the Bose Hubbard and paired-Fermi Hubbard lim- 
its are obtained naturally from the ansatz of Eq. ^ in 
the limits of large negative and positive detuning. The 
FBHH we have proposed therefore correctly obtains the 
endpoints of the BCS-BEC crossover on a lattice. 

Next consider the case of the physically reasonable pa- 
rameter set Vb = Vf, Jb = J'f, with the scaling chosen 
such that Vf ~ 1. The quartic equation has four roots. 
Two are complex and therefore physically extraneous. 
The other two represent an energy minimum and an en- 
ergy maximum. There is no saddle point. The phase 
diagram is shown in Fig. |2]for v = —10,-1, 1/2, 10 and 
g ^ 1. The results are qualitatively the same for all 
g ^ 0- The point = 1/2 is the actual crossover in our 
model, i.e., the point at which the Mott borders become 
degenerate. To illustrate this, in Fig. Ola) is shown the 
mixing angle x as a function of ly. Note the appropriate 
1/ ±oo limits. In Fig. Efb) are shown the y-intercepts 
of the Mott phases from the phase diagrams of Fig. [3 as 
a function of v. These go through an avoided crossing at 
1/ ~ 1/2. Smaller values of \g\ cause the avoided cross- 
ing to become narrower. Similarly, the width of x(^) in 
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FIG. 3: (color online) (a) The mixing angle x as a function 
of the detuning u/Vf. (b) The j/- intercepts in the phase dia- 
gram of Fig.|21go through an avoided crossing as a function of 
the detuning. Blue solid curves: energy minima; red dashed 
curves: energy maxima. 



Fig- EI 8') is proportional to l^j. 

In conclusion, we have proposed a general Fermi-Bose 
Hubbard model which describes the BCS-BEC crossover 
on a lattice. Our restriction of the Hilbert space to the 
lowest band and paired fermions corresponds to the ex- 
perimentally realizable case of from zero to two fermions 
per site in three dimensions and a strongly confining lat- 
tice. We used a superposition ansatz (Eq. which is 
relevant to both broad and narrow Feshbach resonances, 
i.e., for general coupling g. We found that the Paired- 
Fermi Hubbard and Bose Hubbard phase diagrams ap- 
pear naturally and asymptotically for large positive and 
negative detuning. We also showed that the Mott phases 
of the dressed fermion and boson fields go through an 
avoided crossing as the system approaches resonance. 
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